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For a nitrogen dimer in insulating MgO, a ferromagnetic coupling between spin-polarized 2p-holes 
is revealed by calculations based on the density functional theory amended by an on-site Coulomb 
interaction and corroborated by the Hubbard model. It is shown that the ferromagnetic coupling is 
facilitated by a T-shaped orbital arrangement of the 2p-holes, which is in its turn controlled by an 
intersite Coulomb interaction due to the directionality of the p-orbitals. We thus conjecture that 
this interaction is an important ingredient of ferromagnetism in band insulators with 2p dopants. 
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Ferromagnetic (FM) insulators offer the potential for 
use as active barrier materials in novel spintronic tun- 
neling devices. In the past years we have witnessed a 
new vista in engineering FM insulators not by doping 
nonmagnetic insulators with transition-metal ions as is 
traditionally accomplished in diluted magnetic semicon- 
ductors [H, m, but by doping with sp-elements. The 
novel magnetic materials design was encouraged by un- 
expected experimental observations of room-temperature 
ferromagnetism in defective wide-gap oxide semiconduc- 
tors and insulators [sl-Q- This phenomenon (referred to 
as sp- or d^-magnetism) was associated with the partially 
filled p-states of the intrinsic defects like cation/anion va- 
cancies [l, [sl, [sl or first-row (2p) dopants 

Shortly after, numerous cation-deficient or N- and C- 
doped oxides and sulfides were theoretically predicted 
to be FM at room temperature [lol-fl^. Among those, 
MgO is distinct as it is certainly the most important bar- 
rier material for magnetic tunnel junctions. Recently, N- 
doped MgO films were experimentally shown to exhibit 
ferromagnetic properties upon thermal annealing [Tij . 

It is commonly concluded that the FM interaction be- 
tween the defects in d^-magnets is due to partially oc- 
cupied spin-polarized defect states, that are sufficiently 
extended to provide an exchange interaction via the 
double-exchange mechanism. Most likely, however, at 
low concentrations the impurity 2p electrons experience 
a strong on-site Coulomb repulsion U because of their 
spatial localization [l6|, leading to insulating behavior 
that changes the mechanism of exchange interaction and 
could weaken or change the sign of magnetic coupling. 

Correlation effects in defect-free oxides with part ially 
filled oxygen p-shells were studied in Refs. |17H19|. 
The importance of electron correlations in the impu- 
rity p-states was recently investigated using density func- 
tional theory (DFT) approaches amended by an on-site 
Coulomb repulsion [ll|, |20| for the example of N-doped 
MgO. Jahn- Teller-like (JT) distortions, captured by 
2l| or self-interaction correction [22| schemes, were found 
to evoke an energy splitting between the occupied and un- 



occupied nitrogen 2p-states, increasing the localization of 
a spin-polarized hole at one of the 2p-orbitals by pushing 
it deeper into the band gap and thus breaking the initial 
symmetry of the electronic state. 

In this Letter we present new insights into the role 
of electron correlations in N-doped MgO by carrying 
out DFT calculations where strong correlations are ac- 
counted for by the GGA (generalized-gradient approx- 
imation approach, that are conceptualized on 
the basis of a minimal Hubbard model. We demonstrate 
that the symmetry breaking and the subsequent metal- 
insulator transition occurs even without JT distortions, 
and is thus electronically driven. More importantly we 
find for the insulating state of N-dimers, which are likely 
to be formed in the annealing process, a FM spin coupling 
and a T-shaped 2p-orbital arrangement (OA). Employing 
a Hubbard model we show that a weak intersite Coulomb 
repulsion, that is well accounted for by the DFT, com- 
bined with the on-site Coulomb interaction is responsi- 
ble for this result. While in periodic strongly correlated 
3(i systems the favored spin- and orbital arrangement is 
usually well explained by the Kugel-Khomskii model [24 1 
which does not resort to the intersite repulsion, similar 
arguments applied to N-N dimers in MgO yield an anti- 
ferromagnetic (AFM) coupling, in contradiction with our 
first-principles results. The nearest-neighbor FM interac- 
tion is not sufficient to explain the long-range FM order, 
but it is a necessary condition. Our aim is to shed some 
light on the physics driving this interaction. 

We utilized the full-potential linearized augmented 
plane wave code FLEUR [25] in our GGA+/7 calculations. 
Relaxations were performed with an on-site Coulomb en- 
ergy [7=4.6 eV and Hund's exchange J=1.2 eV applied 
on the 2p-states of O and N. The total energies were com- 
pared [26| for a set of U -values that ranges between 3 and 
6 eV, reported for O 2p-states in transition-metal oxides 
from photoemission and Auger experiments [27| , and for 
two values of J, 0.6 and 1.2 eV (from Hund's exchange for 
atomic N and O j26|, we expect the physically relevant 
J- values in a solid to fall within this range). Supercells 
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of 64 atoms of host MgO, with a Brihouin-zone 2x2x2 
k-mesh are used to study a single N impurity and N-N 
dimers at O sites. N-O-N "dimers" are treated in a 144- 
atomic superceh and sampled at the F-point. Different 
OAs are introduced by initiating a specific occupation of 
the p-orbitals of N in the GGA+/7 density- matrix (28| . 
Calculations are performed for all possible types of OA 
at FM and AFM spin alignment (12 different spin-orbital 
configurations in total) for N-N and N-O-N dimers. 

We first consider the case of a single N impurity substi- 
tuting O in MgO to check whether the symmetry break- 
ing and the metal-insulator transition upon applying a 
finite U is driven electronically or by local lattice distor- 
tions. Figure [TJa) shows the N density of states (DOS) of 
the structurally unrelaxed supercell. Within the GGA, 
N introduces two triply degenerate states, an occupied 
spin-up p-state in the valence band and a partially unoc- 
cupied spin-down state at the Fermi level, E^. A single 
spin-polarized hole is evenly distributed among the three 
N p-orbitals. The situation changes drastically upon ap- 
plying U : the minority 2p-state splits into a doubly de- 
generate occupied and a non- degenerate unoccupied level 
(Figdl^b)). This symmetry breaking occurs even without 
lattice distortions, which clearly demonstrates that it is 
an electronically driven effect. Similar behavior was no- 
ticed for the cation vacancy in ZnO 1291. This is different 
from the previous reports on MgO:N [20|, where the main 
effect of both -\-U or self- interaction corrections resulted 
in symmetry breaking via a JT distortion. 

Next, we examine the relative spin and p-orbital ar- 
rangement for two substitutional N impurities in MgO. 
Recent LSDA studies [30| have shown that N-N impurity 
pairing in a structurally relaxed MgO supercell leads to a 
non-magnetic insulating ground state with a fully occu- 
pied bonding ppa (as well as ppn) and completely empty 
antibonding ppa^ states. However, we show that this is 
no longer the case in the presence of strong correlations. 
We consider the following relative positions of the two N 
atoms: (i) at the nearest neighbor (n.n.) sites, at a dis- 
tance (iN-N = 2.97 A - a N-N dimer; (ii) connected via an 
O atom ((^N-N = 5.83 A ) - a N-O-N dimer, and model, 
correspondingly, the 1st and 2nd n.n. interaction between 
N atoms in the oxygen sublattice of MgO. We omit the 
N-Mg-N configuration from the analysis because we find 
the magnetic interactions to be negligibly small (3l| . 

We define a coordinate system with the Px and Py or- 
bitals pointing toward neighboring N (or O) atoms in 
the xy plane (FigHJc)), while the Pz orbital points out of 
plane, toward Mg. In the considered cases of N-N and N- 
0-N dimers, both N atoms are situated along the x axis, 
thus the largest hopping is of ppa type between neigh- 
boring orbitals, while Py andp^ orbitals are essentially 
non-bonding. The sign and strength of the magnetic cou- 
pling between the two spin-polarized holes is found by the 
total energy difference A£^afm-fm between the AFM 
and FM state, calculated within GGA-\-U for eachp-hole 
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Figure 1: (color online) Spin- and orbital-resolved DOS for a 
single N in MgO calculated in GGA (a) and GGA+[/ (b). (c) 
Minority p-holes spin-density distribution in xy-p\a,ne around 
the N-N dimer, with (xx), T-shaped (xy), and (yy) OA. Blue, 
red, and grey spheres denote Mg, O, and N, respectively, 
(d) Spin- and orbital-resolved DOS of N-atom 1 (2) (solid 
(dashed) lines) of the N-N dimer for the depicted OAs, with 
FM spin alignment; red, blue, green line: p^-, Py-, p^-orbitals. 



OA. To be able to disentangle the electronic mechanism 
of OA from the JT one, for the N-N dimer we performed 
calculations with and without structural relaxations. 

The results are summarized in Table HI For each OA 
the spin configuration (FM or AFM) with the lowest en- 
ergy is shown. Clearly, there are substantial energy dif- 
ferences between different types of OA. Comparing the N- 
N and N-O-N dimers, the rapid decrease of A£^afm-fm 
with the N-N distance confirms previous studies [3l| • Im- 
portantly, for both dimers the FM T-shaped OA {{xy) or 
{xz)) is of lowest energy, moreover for any choice of U 
and J parameters (see Supplemental Material (26|). 

It was suggested [III, [2o| that the introduction of U 
would reduce the magnetic coupling in (i^-magnets. How- 
ever, a precise physical value of the parameter U is not 
known. To check the persistence of the magnetic cou- 
pling, we repeated the calculations for the N-N dimer 
varying U in the range between 0-6 eV (neglecting atomic 
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Table I: The GGA+t/ total energies AE relative to the lowest- 
energy state and exchange interaction A^afm-fm both in 
meV/supercell, calculated for different OA for N-N and N- 
0-N dimers in MgO. Positive (negative) A£^afm-fm corre- 
sponds to FM (AFM) couphng. 
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relaxations), with J/U=0.26. At U = J = 0, the system 
is half- metallic and ferromagnetic with A£^afm-fm = 
265 meV. However, at /7 > 2 eV, it becomes insulating 
and the interaction through the double-exchange mecha- 
nism is no longer possible. Thus, the magnetic coupling 
is much weaker and we expect AFM coupling due to ki- 
netic exchange interaction. Instead, orbital arrangement 
sets in and despite the insulating state, the interaction 
between the spin- 1/2 holes remains FM with total mag- 
netic spin moment of 2/iB- In the following we provide a 
qualitative understanding of this unexpected finding. 

In the first step we show that kinetic exchange alone 
favors an AFM ground state. We consider a minimal 
multiorbital Hubbard model that describes the system of 
strongly correlated open-shell 2]9-electrons in terms of the 
kinetic energy, t, of spin-conserving hoppings between or- 
bitals of largest overlap (n.n. ppa-orbitals, see FiglSfa)), 
the on-site Coulomb repulsion energy U and the interor- 
bital exchange (Hund's rule) coupling J (see also Eq. ([T])). 

For a given OA, the favored magnetic interaction is 
determined by hopping of the electrons within the N- 
N-dimer. A substantial magnetic coupling occurs only 
for those OAs that allow hopping between the two Px- 
orbitals of the N-dimer, where at least one of these or- 
bitals is half-filled: the {xx) and the T-shaped {xy) and 
{xz) OAs (see Table H]). The relative coupling strength 
can be qualitatively understood by considering that each 
spin-conserving hopping t into half-filled orbitals lowers 
the energy by t^/{U — J) (in 2nd order perturbation 
theory) if the intermediate virtual atomic state satisfies 
Hund's rule for one of the atoms, and /U otherwise. 
In the {xx) OA, this gives Eafm = —2t^/U^ Eym = 0, 
i.e. an AFM ground state (the factor 2 arises because 
there is one half-filled orbital per atom participating in 
the hopping). In the T-shaped OA there is only one half- 
filled orbital participating; then we have £^afm = —t^/U^ 
Eym = —t^/{U — J), i.e., a FM ground state. 

Among all possible types of OA, the T-shaped and 
AFM {xx) states are the ones that have the lowest en- 
ergy, because these are the only ones that allow hopping 
between the N atoms. The hopping between O and N 



atoms provides a smaller total energy gain, because of (i) 
the separation between the on-site energy levels of O and 
N, that differ by 0.5-1 eV, and (ii) the delocalization of 
the O states that form an itinerant band. From the previ- 
ous discussion follows that among the two OAs, the AFM 
{xx) one should be lower in energy by 2t^ /U — /{U — J) 
as long as J < U/2. Thus, the arguments of the kinetic 
exchange interaction favor the AFM {xx) state, contra- 
dicting our DFT finding of a FM T-shaped ground state. 

Obviously, a qualitative description of the relative sta- 
bility of spin- and orbitally-arranged states in N-doped 
MgO requests an extension of the minimal Hubbard 
model beyond the conventional electron hopping t and 
on-site Coulomb repulsion U. The N-N dimer breaks 
the cubic symmetry of the MgO lattice. Considering 
the directional nature of the 2p-orbitals, this leads to 
an orbital-dependent intersite Coulomb repulsion, that 
is well accounted for within the GGA/LDA. In terms of 
the Hubbard model it is expressed by an additional effec- 
tive intersite Coulomb repulsion energy V between n.n. 
ppcr-orbitals of N-N and N-0 pairs. To capture this effect 
by such an extended Hubbard model we study an impu- 
rity cluster of 2 N and 6 n.n. O atoms (Figl2fa)), with 
1 (2) p-orbital on each O (N) site that form n.n. ppcr- 
orbitals. Thus, we restrict ourselves to 10 p-orbitals that 
can be occupied by 18 electrons and 2 holes of spins t and 
I. The Hamiltonian for such a system can be written as: 

771,3 <7n,7n'> s 

+ ^ ^ n^t^mi Y nmUm' (1) 

771 <m,m'> 

H~ {U J) ^ ^ Tims'^m' s H~ U ^ ^ Tlms'^m' — s 
^m,m'^,s <^m,77i''^,s 

Here <m, m^> indicates the pair of n.n. orbitals of ppa- 
overlap, Cm denotes the energy level, cj^g {cms) creates 
(annihilates) a particle of spin s =t, | of orbital m, while 
^ms, = ^mt + ^mi are the spin- number and num- 
ber operators, respectively. Since we have two orbitals 
at the N atoms, accounts for the mter-orbital on-site 
Coulomb interaction; <^m^m'^ stands for the two on- 
site orbital pairs (1,2) and (3,4). We neglect on-site spin- 
flip and pair-hopping terms as they are non-zero only for 
the states with an empty and a fully occupied orbital at 
the same site, which turned out to be of much higher 
energy. The matrix elements of ([T]) are evaluated in the 
18-electron-state basis (190 basis functions in total). The 
eigenvalues of the model are obtained numerically by ex- 
actly diagonalizing ([1]) for a given set of parameters. We 
impose the relations U > U' > J>0, = U - 2J 
(for a derivation see the Supplemental Material 0]) and 
we take t = 0.75 eV (extracted from the splitting of the 
DOS on FiglU^d)). As the p- levels of N atoms are higher 
in energy than those of O, we take = 0.5 eV for N and 
zero for O. The p^-orbitals which are omitted from the 



4 



model are considered to be always filled and have only 
an effect of a constant shift on the on-site energies e^. 
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Figure 2: (color online) (a) A schematic view of the clus- 
ter of p-orbitals, used in the model Hamiltonian ([T|). The 
dashed blue line depicts px- and p^-orbitals of the n.n. nitro- 
gen atoms (Ni and N2), while the solid black line depicts the 
corresponding orbit als of the O atoms. The possible hopping 
processes t between two orbitals on neighboring sites are also 
shown, (b) The {U — 2 J) vs. U phase diagram of the model 
for V = and V = 0.2 eV. Orange circles denote the U and 
J values, considered in our ab initio calculations [2^, outlin- 
ing the physically relevant area (in green), (c) The phases I, 
II, III and IV, corresponding to different sequences of several 
orbitally-arranged states in the order of energetic stability. 



Figure [2fb) shows a {U — 2J) vs. U phase diagram, re- 
sulting from Hamiltonian ([T]) for V = and V = 0.2 eV. 
The black lines denote the boundaries between the phases 
I, II, III and IV, corresponding to different sequences 
of energetic stability of OAs (see Fig. El^c)). We find 
that the FM (xy) state is the lowest-energy one for the 
phases II, III and IV. However, for almost all U and J 
values, the AFM (xx) state (phase I) is the ground state 
if F = 0, contradicting our ab initio results [26|. More- 
over, at V = the sequence of DFT states from Table |I] 
(phase HI) is not reproduced. In contrast, at finite V the 
whole sequence of DFT states is reproduced for a wide 
range of physically relevant values of U and J. Hence, 
the inter site Coulomb interaction V, which defines the 
repulsion of electrons on n.n. p-orbitals directed towards 
each other, is crucial for stabilizing the FM state in doped 
insulators with partially occupied p-orbitals like MgO:N. 

In summary, we have shown that the physics of spin- 
polarized holes of N-doped MgO is largely determined by 
p-electron correlation effects. For a single impurity, the 
splitting of the N 2p-states at finite U is mainly an elec- 
tronic effect, enhanced by the lattice distortions. More- 
over, the p-electron correlations lead to a spin and or- 
bital arrangement of the 2p-hole states of the N-N dimer, 
resulting in a FM coupling between T-shaped orbitally- 
arranged spin-polarized p- holes. The FM state is realized 
at reasonable values of U and J, if an intersite Coulomb 
repulsion V is accounted for. Considering the small value 



oi V ^ 0.2 eV, necessary to drive the phase transition 
from an AFM to a FM ground state, and the fact that the 
directional nature of the p-orbitals controls the strength 
of the intersite Coulomb repulsion, we conjecture that 
our finding affects all oxides where magnetism is due to 
p-electrons or defects in the p-electron systems. 
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